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We perform a global analysis of all recent experimental data from elastic parity- violating electron 
scattering at low Q 2 . The values of the electric and magnetic strange form factors of the nucleon 
are determined at Q 2 = 0.1 GeV/c 2 to be G% = -0.008 ± 0.016 and G S M = 0.29 ± 0.21. 
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The existence of a "sea" of quarks and antiquarks in 
the nucleon has been firmly established in deep-inelastic 
lepton scattering experiments as well as in the produc- 
tion of dilcpton pairs (the Drell-Yan process). However, 
demonstrating the role of these qq pairs in the static elec- 
tromagnetic properties of the nucleon has been a more 
elusive and difficult task. 

As the lightest quark that contributes only to the qq 
sea, the strange quark provides a unique window on the 
role of the sea in the nucleon's electromagnetic structure. 
As suggested by Kaplan and Manohar [1|, knowledge of 
neutral current form factors, when combined with elec- 
tromagnetic form factors, provides access to the contri- 
bution of strange quarks to these form factors. At low 
momentum transfers, the neutral current form factors 
can be determined through parity- violating (PV) elec- 
tron scattering experiments [2, El • 

During the last decade, there has been dramatic 
progress in the study of the strange quark-antiquark con- 
tributions to the nucleon elastic electromagnetic form 
factors. A series of definitive PV electron scattering ex- 
periments along with several theoretical studies now pro- 
vide a basis for extracting precision information on these 
strange quark contributions. In this work we report the 
results of a global analysis of all these experiments, in- 
cluding both the latest data obtained in experiments per- 
formed at the Jefferson Laboratory and appropriate the- 
oretical input on radiative corrections, and obtain values 
for the strange electric and magnetic form factors of the 
nucleon at a four-momentum transfer Q 2 = 0.1 GeV/c 2 . 
We have also studied the sub-leading Q 2 dependence of 
these two form factors, and find that so far the data do 
not provide conclusive information. 



I. STRANGE FORM FACTORS AND 
PARITY- VIOLATING ELECTRON SCATTERING 

The nucleon vector strange form factors, G E and G S M , 
characterize the contribution of the strange sea quarks 
to the nucleon electromagnetic form factors, and thereby 
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their contribution to the charge and magnetization dis- 
tributions in the nucleon. With polarized electron fa- 
cilities, G S E and G S M can be accessed by measuring the 
PV asymmetries in elastic e-p scattering ,_quasielastic e-d 
scattering, and elastic e- 4 He scattering In very gen- 
eral terms, the parity- violating asymmetry Apy can be 
written as 



A PV = A nvs + rj E G s E + r] M G 



M • 



(1) 



where A nvs is the "non- vector-strange" asymmetry (inde- 
pendent of G S E and G S M ), and r] E and t]m are functions of 
kinematic quantities, nucleon electromagnetic form fac- 
tors, and nuclear models (for non- hydrogen targets). 

For elastic e— p scattering, the full form of the asym- 
metry is [H 



G F Q 2 



1 
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x {(e(G E f + r(G p M ) 2 )(l - 4sin 2 W )(1 + R p v ) 
- (eG p E G E + tG p m Gm){1 + R v ) 



(eG p E G s E + TG p M G s M )(l + R^) 
e'(l - 4sm 2 w )G p M G e A } , 



(2) 



with 



AM 2 



e= 1 + 2(1 + T)tan 2 - 



e'= Vr(l+r)(l-e2), 

where M p is the mass of the proton and 9 is the electron 
scattering angle. In Eqn. [2j Gf and a are the Fermi 
and fine structure constants, respectively. Q 2 is the four 
momentum transfer. G E '^ are the proton and neutron 
electric and magnetic form factors, while G e A is proton 
axial form factor seen by an electron. In order to extract 
contributions from G S E M to A P PV , one must include the 
effects of Standard Model (SM) 0(a) electroweak radia- 
tive corrections @ ■ It is often useful to characterize these 
corrections in terms of ratios i?y",A of the 0(a) hadronic 
vector (V) and axial vector (A) weak neutral current am- 
plitudes to the corresponding tree-level amplitudes. The 
Ry, R v , Ry^ give these ratios for vector proton, neutron, 
and SU(3)-singlct amplitudes, respectively. In principle, 
their values can be obtained using the SM predictions 
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for the effective electron-quark couplings C\ q given in 
Q. However, the quoted C\ q do not include perturbative 
QCD contributions or coherent strong interaction effects 
in the radiative corrections associated with elastic scat- 
tering from a nucleon. A recent analysis of these effects 
has been given in Ref. Q and up-dated in Ref . ■ The 
latter work also gives an improved treatment of strong in- 
teraction contributions to the running of the weak mixing 
angle in the M S renormalization scheme from its value at 
the Z-polc, s 2 z = sin 2 6\y(Mz), to the quantity appropri- 
ate for precise, low-energy neutral current experiments, 
sin 2 ^iy(0). All of these effects are included in the Ry 
given in Table HI The theoretical uncertainties in Ry 

and Rp are less than one percent and have a negligible 
impact on our analysis, so we do not quote these errors in 
Table |U The theoretical error in (1 - 4 sin 2 6 W ){1 + R P V ) 
is ±0.0008 [Tc|, or slightly more than one percent. Since 
this error receives roughly equal contributions from the 
uncertainty in s 2 z and from the 0{a) Z7 box graph cor- 
rections, it is not appropriate to quote an uncertainty in 
Ry alone. The uncertainties associated with other strong 
interaction effects are sub-dominant. For the kinematic 
range of the PV experiments analyzed here, the Ry have 
a negligible impact on the Q 2 -dependence of A p pv and 
are taken to be constant. We use the conventional MS 
renormalization scheme. Therefore, sin 2 9w m Eqn. [2] 
and hereafter shall take its value as s 2 z . 

For reader's convenience, constant parameters (in some 
cases also the uncertainties) appearing in Eqn.[2]and later 
in Eqns.[3]and[2]are summarized in Table [IJ 



Parameter 


Value 


a 


1./137.03599911 


4 


0.23122(15) 


G F 


1.16637 x 10- 5 /GeV 2 


M p 


0.98272 GeV 


Ai 


1.00(0.04) (GeV/c) 2 



Parameter Value 



9a/ gv 
3F-D 
As 

Ry 
I?" 



R 



(0) 



-1.2695 
0.58(0.12) 
-0.07(0.06) 
-0.0520 
-0.0123 
-0.0123 



TABLE I: Summary of parameters (with some uncertainties 
in parentheses) in Eqns. [2] [3] and [4] The values of a, s%, 
Gf, M p , qa I gv are taken directly from Q. Ry, and Ry' 1 are 
converted from C\ q parameters in 
proton weak charge given in Ref. [1 
and uncertainty of from [TJ], 3-F 
from [13]. 



(0) 
V 

Ry is derived from the 
. We adopt the value 
- D from [131. and As 



The effective axial form factor G\ receives a number 



of contributions and may be written as 



G A (Q 2 ) = G D (Q 2 



£d ( i 

l9v 



where 



R 



T=l\ 



IX 

3F 



D 



-R 



T=0 



A s (l + 4 0) ) 



G D (Q 2 
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QV^A? 



(3) 



(4) 



parameterizes the <5 2 -depcndence with a dipole form and 
h-A is the corresponding axial dipole mass. The ratio 

— — is the isovcctor axial form factor of the nucleon at 

av 

zero momentum transfer, which is precisely measured in 
the neutron beta decay. F and D are the octet baryon 
beta decay parameters, which can be determined by com- 
bining data from neutron and hypcron beta decays un- 
der the assumption of SU(3) flavor symmetry. As is 
the strange quark contribution to nucleon spin. Assum- 
ing a gentle evolution from the perturbative to the non- 
perturbative domain, this quantity can be obtained from 
inclusive, polarized deep inelastic lepton-nucleon scatter- 
ing. 



The ratios 



R T A=° 



', and R A characterize the ef- 
fects of electroweak radiative corrections to the isovector, 
isoscalar, and SU(3) singlet hadronic axial vector ampli- 
tudes. Note that while R A =1 and R A ^ give the ratios of 
the 0(a) and tree-level hadronic axial vector neutral cur- 
rent amplitudes in the isovector and SU(3) singlet chan- 
nels, respectively, R^ =0 does not have this interpretation 
since the tree-level isoscalar hadronic axial vector ampli- 
tude vanishes in the SM. 

Conventionally, these quantities are divided into two 
pieces: the "one-quark" and "many-quark" contribu- 
tions. The one-quark contributions correspond to renor- 
malization of the effective vector electron-axial vector 
quark couplings, Ci q , and their values can be obtained 
from the SM predictions for these couplings given in 
Ref. @. The many-quark contributions include the so- 
called "anapole" effects as well as coherent strong inter- 
action contributions to the radiative corrections. In con- 
trast to the situation with the vector corrections, Ry , the 
relative importance of many-quark effects in the Ra can 
be quite pronounced. The small vector coupling of the 
electron to the Z-boson, gy = —1 + 4 sin 2 §w ~ —0.075, 
leads to a suppression of the tree-level hadronic axial vec- 
tor amplitude. However, gy is absent from a variety of 
both one- and many-quark radiative corrections. Thus, 
one would expect the magnitudes of the Ra to be of or- 
der several percent, rather than the generic a/w ~ 0.3% 
scale normally associated with electroweak radiative cor- 
rections. As a result, the impact of otherwise negligible 
strong interaction effects in the many-quark corrections, 
such as the anapole contributions, can be amplified [3]. 

An appropriate framework for treating the many-quark 
effects associated with physics at the hadronic scale is 
chiral perturbation theory. A comprehensive analysis of 



and -R^ -0 nas been 



15| . This analysis 



the anapole contributions to R A =1 
carried out to chiral order p 3 in Ref. 
included both one-loop contributions associated with the 
octet of pseudoscalar mesons as well as the full set of low- 
energy constants (LECs) that arise at this chiral order. 
A generous theoretical range was assigned to the LECs, 
leading to a quoted theoretical uncertainty in the total 
Ra that is larger than the (logarithmically enhanced) 
one-quark corrections. The theoretical SM uncertainty 
is likely to be smaller. 
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The corresponding results, updated for the present 
value of the weak mixing angle, are given in Table [U The 
resulting prediction for G e A is consistent with both the 
results of the SAMPLE deuterium measurement [lj| [TtJ , 
which is particularly sensitive to the dominant isovector 
axial component, as well as other theoretical models for 
the anapole contributions [HI, . No evaluation of the 
"many-quark" contribution of R^ A has been made in the 
literature. We assume it is zero and assign the size of 
"one-quark" value for R ?' as its uncertainty. 



the first and second uncertainties for the G° data are the 
uncorrelated and correlated experimental uncertainties, 
respectively. The values of r] E and t\m are also listed 
in the table. In calculating them, we have adopted a re- 
cent parametrization of the nucleoli electromagnetic form 
factors from Ref . . For the SAMPLE deuterium mea- 
surement, the t\m is taken from Ref. [l6j |. whereas its r\ E 
is taken to be 1.79 according to the static approximation. 





Ra 1 


H A 


D 

K A 


one-quark 


-0.172 


-0.253 


-0.551 


many-quark 


-0.086(0.34) 


0.014(0.19) 


N/A 


total 


-0.258(0.34) 


-0.239(0.20) 


-0.55(0.55) 



III. 



GLOBAL ANALYSIS 



TABLE II: The "one-quark" and "many-quark" [H cor- 
rections to the axial charges, both in MS, as well as the com- 
bined corrections. 

The PV asymmetry for the neutron can be obtained 
by exchanging the "p" and "n" indices on nucleon form 
factors in Eqn.[51 and flipping the sign of the first isovec- 
tor term in the expression for G\ in Eqn. [5] To first 
order, the PV asymmetry from a deuterium target is a 
cross-section weighted average of the proton and neutron 
asymmetries, which leads to an enhancement of the con- 
tribution of G\ and a suppression to the relative contri- 
bution due to G S E and G S M . Obviously a nuclear correc- 
tion needs to be applied in the analysis. In this note, for 
the SAMPLE deuterium measurement, we shall adopt 
the asymmetry expression given in . 

The He nucleus is spin zero, parity even and isoscalar. 
The P V asymmetry takes a much simpler form Q : 



A 



He 
PV 



G F Q 2 

(4 sin 2 6W(l + i?y = °) + 



2(1 + R i ° ) )G 



, (5) 



where the isoscalar Ry factor is related to RS, and R\\ 



B T=0 



R v - (1-4 sin 2 e w )R v 



4 sin 2 9 



(0) 



ti- 



ll. 



EXPERIMENTAL DATA 



In this section, the world data of PV clastic scatter- 
ing within a Q 2 range from 0.07 to 0.5 (GeV/c) 2 will 
be summarized. These include SAMPLE-H [13, 
SAMPLE-D [H 03], HAPPEx-H-99 [HI], HAPPEx-H- 
a [H, HAPPEx-He-a [2l, HAPPEx-H-b and HAPPEx- 
He-b 0, PVA4-H-a 0, PVA4-H-b [H[, and the first 
14 Q 2 bins in G° forward angle [27||. The kinematics, 
targets, and the measured asymmetries in these exper- 
iments are summarized in Table IIII1 In column A p h ys , 



A. A nV s and Theoretical Uncertainties 

We shall now present a combined analysis of the world 
data aiming to extract G S E and G S M . A global fit, gen- 
erally speaking, is obtained by simultaneously solving a 
set of equations 



m l ± <x(m,) = t i (a 1 ,a,2, - ■ •) ± cr(ti) 



(7) 



where and U(ai, a%, ■ ■ • ), respectively, are the mea- 
sured and theoretical values for experiment i. In this 
expression, a{mi) and a(ti) are their uncertainties, and 
01, 02, • • ■ are the free parameters one seeks to determine. 
In our case, 



A 1 + 



A n 



i] E G s E + r) M G 



M ) 



(8) 



with G S E and 67^ being the free parameters. In the pre- 
vious section, we have discussed the value and uncer- 
tainty of A l phys , as well as t\e and r\M (Table ITll) . For 
each measurement, the values of A nvs can be also com- 
puted straightforwardly using the formalism in Sec. |U 
They are listed in Table |lVl We again have used the 
parametrization of the nucleon electromagnetic form fac- 
tors from Ref. [29[. As mentioned, the A nvs for the SAM- 
PLE deuterium measurement is calculated based on the 
asymmetry expression in fiH ] with the theoretical value 
ofG%. 

The treatment of the theoretical uncertainties a(ti) is 
more subtle. a(ti) receives dominant contributions from 
the following sources: the nucleon axial form factor {G e A ), 

nucleon electromagnetic form factors {G E '^J), and nu- 
clear corrections. Theoretical uncertainties from a given 
source are correlated among different experiments. The 
uncertainty in G e A can be calculated based on the un- 
certainties in Tables U and [H and is dominated by the 
uncertainty of the "many-quark" electroweak radiative 
corrections on the Ra factors in Table |TJ For the nucleon 
electromagnetic form factors, based on the spread of the 
world data (see, e.g. Ref. [29[), we estimate their relative 
uncertainties as 



u E 



2.5% 



CP 



= 1.5% , 



= 15% and 



°"(Gm) 

fin 



= 1%, 
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Experiment 


Target 


Q 2 

(GeV/c) 2 


Olab 

(deg) 


(ppm) 


'IE 

(ppm) 


(ppm) 


Reference 


SAMPLE-H 


H 2 


0.1 


144.4 


-5.61 ± 1.11 


2.07 


3.47 


[17, 20] 


SAMPLE-D 


D 2 


0.091 


140.8 


-7.77 ± 1.03 


1.79 


0.77 


[16, 12] 


HAPPEx-H-99 


H 2 


0.477 


12.3 


-15.05 ±1.13 


56.89 


22.62 


[21] 


HAPPEx-H-a 


H 2 


0.099 


6.0 


-1.14 ±0.25 


9.55 


0.76 


[22] 


HAPPEx-H-b 


H 2 


0.109 


6.0 


-1.58 ±0.13 


10.59 


0.93 


[24] 


HAPPEx-He-a 


4 He 


0.091 


5.7 


6.72 ±0.87 


20.19 





[23J 


HAPPEx-He-b 


4 He 


0.077 


5.8 


6.40 ± 0.26 


16.56 





[2J 


PVA4-H-b 


H 2 


0.108 


35.52 


-1.36 ±0.32 


10.08 


1.05 


[26] 


PVA4-H-a 


H 2 


0.23 


35.45 


-5.44 ±0.60 


22.56 


5.07 


[25] 


G° 


H 2 


0.122 


6.68 


-1.51 ±0.49 ±0.18 


11.96 


1.17 


[27] 


G° 


H 2 


0.128 


6.84 


-0.97 ±0.46 ±0.17 


12.60 


1.30 


[27] 


G° 


Ho 


0.136 


7.06 


-1.30 ±0.45 ±0.17 


13.46 


1.48 


[27] 


G° 


Ho 


0.144 


7.27 


-2.71 ±0.47 ±0.18 


14.32 


1.66 


[27] 


G° 


Ho 


0.153 


7.50 


-2.22 ±0.51 ±0.21 


15.31 


1.89 


[27] 


G° 


Ho 


0.164 


7.77 


-2.88 ±0.54 ±0.23 


16.53 


2.19 


[27] 


G° 


H 2 


0.177 


8.09 


-3.95 ±0.50 ±0.20 


17.99 


2.58 


[27] 


G° 


Ho 


0.192 


8.43 


-3.85 ±0.53 ±0.19 


19.69 


3.07 


[27] 


G° 


Ho 


0.210 


8.84 


-4.68 ±0.54 ±0.21 


21.77 


3.71 


[27] 


G° 


Ho 


0.232 


9.26 


-5.27 ±0.59 ±0.23 


24.37 


4.60 


[27] 


G° 


Ho 


0.262 


9.92 


-5.26 ±0.53 ±0.17 


28.00 


5.99 


[27] 


G° 


Ho 


0.299 


10.63 


-7.72 ±0.80 ±0.35 


32.60 


7.99 


[27] 


G° 


Ho 


0.344 


11.46 


-8.40 ±1.09 ±0.52 


38.40 


10.89 


[27] 


G° 


H 2 


0.410 


12.59 


-10.25 ±1.11 ±0.55 


47.28 


16.10 


[27] 



TABLE III: A summary of the world data on PV elastic electron scattering within the range of 0.07 (GeV/c) 2 < 
Q 2 <0.5 (GeV/c) 2 , including the average kinematics, targets, published asymmetries A p h yB , as well as coefficients r/E and 
r]M. A p hys, We and T)m are in units of parts-per-million (ppm). The central kinematics for the two PVA4 measurements are ob- 
tained from [28| • For Aphys , the first and second uncertainties for the G° data are the uncorrelated and correlated experimental 
uncertainties, respectively. 



respectively. This is consistent with the uncertainty as- 
signment in Ref. [22|, except that the uncertainty of G E 
in [22j was assigned more conservatively to be 30%. Also 
note that we have made a simplifying assumption that 
these form factor uncertainties are "scaling" in nature, 
independent of the Q 2 . For an analysis with relatively 
small range of Q 2 , this is reasonable. Nuclear correc- 
tions are only relevant for non-hydrogen targets. For the 
two 4 He measurements, according to Ref. [24J, 3% is as- 
signed as the fractional theoretical uncertainty of A nvs . 
Nuclear corrections for the SAMPLE deuterium measure- 
ment have little impact on the final fit, and are there- 
fore neglected. The theoretical uncertainties for the U in 
Eqn.[S]duc to the different sources are summarized in the 
last six columns in Table ITVl To be precise, the content of 
each column gives the change in A nvs ± rj E G s E ± timG s m 

when the source magnitude (|G^|, \G E '^\, or nuclear 
correction to 4 He data) is increased by one standard devi- 
ation. Notice that the nucleon electromagnetic form fac- 
tors also affect the value of rj E and tjm, therefore generate 
"pull terms" linear in G S E and G S M in Table [TV] Such pull 
terms are neglected for the SAMPLE deuterium measure- 
ment. 



B. Combined Analysis at Q 2 =0.1 (GeV/c) 2 

As seen from Table IIII1 a wealth of data exist with 
Q 2 in the vicinity of 0.1 (GeV/c) 2 , including SAMPLE- 
H, SAMPLE-D, HAPPEx-H-a, HAPPEx-H-b, HAPPEx- 
He-a, HAPPEx-He-b, PVA4-H-b, and low Q 2 data from 
G°. It is natural to first make a combined analysis at 
Q 2 = 0.1. To interpolate all data to a common Q 2 , 
we assume G S E oc Q 2 and G S M is a constant [3(j. That 
is, we replace veG s e (Q 2 ) + r? M G s M (Q 2 ) by fj B G E {0.1) ± 

VmG s m (0A), where rf E = 7] E 01 ( q cV/c)2 ■ To simplify our 
notation, we shall use G S E and G S M hereafter to denote 
their values at Q 2 =0A (GeV/c) 2 . In the {G S E , G S M ) space, 
each measurement i provides a linear constraint as 

VeG s e + VmG 6 m = A l phys — A l nvs . (9) 

In Fig. [T] each constraint is shown as a linear band 
in the (G S E ,G S M ) plane, where <?(A phys ) (see Table HTl|) 
and the theoretical uncertainty a(ti) (see Table ITV| have 
been combined in quadrature into an overall uncertainty. 
Somewhat arbitrarily, we include the 3 lowest Q 2 bins 
from G° data in this part of the analysis. For visual 
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Exp 


Q 2 












Theoretical uncertainty 


t(U) (ppm) 








(GeV/c) 2 


ppm 


G\ 




M 




^ M 


N.C. 


SAMPLE-H 


0.1 


-6.85 


0.57 


0.06 + 0.03G B 




0.03G„ 


0.06 - 


- 0.05Gb 


- 0.03Gm 


0.01 


— 0.05 





SAMPLE-D 


0.091 


-8.37 


0.60 




0.05 








0.02 




0.02 


— 0.03 





HAPPEx-H-99 


0.477 


-15.96 


0.23 


0.38 + 0.07GI; 




0.52Gm 


-0.03 


- 0.90Gb 


- 0.02Gm 


0.49 


-0.16 





HAPPEx-H-a 


0.099 


-1.42 


0.01 


0.04 


-0.15GJ, 




0.03Gm 


-0.01 


- 0.05Gb 


+ 0.01Gif 


0.05 


-0.01 





HAPPEx-H-b 


0.109 


-1.64 


0.02 


0.05 


-0.16G B 




0.04G1, 


-0.01 


- 0.06Gb 


+ 0.01G|f 


0.06 


-0.01 





HAPPEx-He-a 


0.091 


7.50 







-0A7G% 









-0.12G|; 





0.23 


HAPPEx-He-b 


0.077 


6.35 







-0.39G| 











-0.08Gb 





0.19 


PVA4-H-b 


0.108 


-2.02 


0.09 


0.06 


-0.14Gb 





0.04GI/ 


-0.02 


- 0.07Gb 


+ O.OIGh 


0.05 


-0.01 





PVA4-H-a 


0.23 


-6.26 


0.26 


0.17 


- 0.13G B 





0.15G^ 


-0.03 


- 0.26G! 


+ 0.02G S M 


0.18 


-0.05 





G° 


0.122 


-1.94 


0.02 


0.06 


- 0.17Gb 


— 


0.04Gm 


-0.02 


- 0.08Gb 


+ 0.01Gif 


0.07 


-0.02 





G° 


0.128 


-2.09 


0.02 


0.06 


-0.17Gb 




0.05Gm 


-0.02 


- 0.09Gb 


+ 0.01Gif 


0.07 


-0.02 





G° 


0.136 


-2.29 


0.02 


0.07 


- 0.18Gb 




0.05Gm 


-0.02 


-0.10Gb 


+ O.OIG^ 


0.08 


-0.02 





G° 


0.144 


-2.50 


0.03 


0.07 


-0.18Gb 




0.06Gm 


-0.02 


-0.11Gb 


+ O.OIG^ 


0.09 


-0.02 





G° 


0.153 


-2.74 


0.03 


0.08 


-0.19Gb 




0.07Gm 


-0.02 


-0.12Gb 


+ 0.01G^ 


0.10 


-0.02 





G° 


0.164 


-3.05 


0.04 


0.09 


-0.19Gb 




0.08Gm 


-0.02 


-0.13Gb 


+ O.OIG^ 


0.11 


-0.03 





G° 


0.177 


-3.43 


0.04 


0.10 


- 0.19Gb 




0.09Gm 


-0.03 


- 0.16Gb 


+ 0.02G1, 


0.13 


-0.03 





G° 


0.192 


-3.88 


0.05 


0.11 


-0.19GI, 




O.IOGm 


-0.03 


-0.18Gb 


+ 0.02GJ, 


0.14 


-0.03 





G° 


0.210 


-4.46 


0.06 


0.13 


- 0.19Gb 




0.12G^ 


-0.03 


-0.21Gb 


+ 0.02G S U 


0.16 


-0.04 





G° 


0.232 


-5.21 


0.07 


0.15 


-0.19Gb 




0.15Glf 


-0.04 


- 0.25Gb 


+ 0.02G S M 


0.19 


-0.05 





G° 


0.262 


-6.30 


0.08 


0.18 


-0.18GI, 




0.18Gm 


-0.04 


- 0.32Gb 


+ 0.02G S M 


0.23 


-0.06 





G° 


0.299 


-7.75 


0.11 


0.21 


-0.15Gb 




0.23Gm 


-0.04 


- 0.40Gb 


+ 0.02G S M 


0.27 


-0.07 





G° 


0.344 


-9.65 


0.14 


0.26 


-0.11Gb 




0.30G S M 


-0.04 


-0.51Gb 


+ 0.02G S M 


0.33 


-0.09 





G° 


0.410 


-12.70 


0.19 


0.32 


- 0.03Gb 




0A0G 3 M 


-0.04 


- 0.69Gb 


+ O.OOGif 


0.41 


-0.12 






TABLE IV: A summary of A nva and the theoretical uncertainties of individual measurements. The uncertainties are grouped 
by six different sources: G\, G P E , G P M , G B , G^ and the nuclear corrections (N.C.) to 4 He data. See text for details. 



clarity, they are combined into a single constraint as 

^ P hys ~ Kvs = 0.84 ± 0.34 = 16.38G| + 1.32G^ , (10) 

which is shown as the solid brown band in Fig. [T] From 
the figure, one sees that the agreement among different 
measurements is generally good. The G° and PVA4 ap- 
pear to be offset from the HAPPEx-H-b measurement, 
but they nevertheless agree within 2a. As explained in 
Sec.|H the SAMPLE deuterium measurement (dashed red 
band) has much less sensitivity to G S E and G S M . 

The 10 different measurements (3 from G°, and the 
other 7 from separate experiments) above provide re- 
dundancy in the joint determination of (G S E ,G S M ). To 
solve for (G S E ,G S M ) and determine the confidence con- 
tours, we follow the standard least square procedure (see, 
e.g. Ref [3l|). Specifically, we rearrange Eqn.[9]into the 
form of Eqn. [7] as 

6 

mi = ti(G%, G S M ) ± a + J2 ± &,k ■ (11) 

fe=0 

where m., and ti are given in Eqn. [51 and Ui is the uncor- 
rected experimental uncertainty. (3i^ denotes the corre- 
lated uncertainty for measurement i with "source index" 
fc. In our case, /J^o equals the correlated experimental 
uncertainty for the G° data and for other experiments, 



and /3i,fc=i.2,3,4,5,6 are the correlated theoretical uncer- 
tainties for each measurement i due to different sources 
(Table [IV|. Then for each given pair of (G S E ,G S M ), the 
X 2 is calculated as 

x 2 = EE( m * - tiW^hi^ - tj) , (12) 

i J 

where i and j arc indices of the measurements, and V 

6 

is the variance matrix with Vij = crfSij + Pi,kPj,k- 

k=Q 

It has been shown in Ref. [31| that the \ 2 constructed 
this way satisfies the standard x 2 distribution, and the 
solution (best fit) can be found by minimizing this x 2 - 
Applying this technique to the 10 measurements in Fig.[TJ 
we obtain 

G% = -0.006 ±0.016, G" M = 0.33 ±0.21 (13) 

with a correlation coefficient of —0.83 between the two, 
and a minimum Xmin = 9-90 with 8 degrees of freedom. 
Note that the uncertainties above are la (Ax 2 = 1) 
"marginalized" uncertainties corresponding to the pro- 
jections of the error ellipse onto the two axes. That is, 
for a given value of G S M , the \ 2 is minimized by vary- 
ing G S E and vice versa. The range defined by this un- 
certainty for a given parameter corresponds precisely to 
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Q"=0.1 (GeV/c) ' 

SAMPLE-H 

SAMPLE-D 

HAPPEx-H-a 

HAPPEx-H-b 

HAPPEx-He-a 

HAPPEx-He-b 

PVA4-H-b 

— G" 




FIG. 1: (Color online) The world data constraints on 
{G%,G S M ) at Q 2 = 0.1 (GeV/c) 2 . The form factors of 
Kelly are used. Different bands in th e p lot represent: 
SAMPLE-H [lljsolid red), SAMPLE-D [Q (dashed red), 
HAPPEx-H-a n23] (dashed blue), HAPPEx-H-b [U (solid 
blue), HAPPEx-He-a [|1 (dashed pink), HAPPEx-He-b 0) 
(solid pink), PVA4-H-b [2Q] (solid green), and the lowest three 
Q 2 bins in G° forward angle [13] (solid brown). The yellow 
and gray blue (dark) ellipses represents 68.27% (Ax 2 = 2.3) 
and 95% (Ax 2 = 5.99) confidence contours around the point 
of maximum likelihood at (G S E = -0.006, G S M = 0.33). The 
black cross represents G S E = G S M = 0. 



to PVA4-H-b), and then systematically included more 
data in the fit with increasing Q 2 . In Table fVl the re- 
sulting (G S E ,G S M ) and the fit quality are summarized. 
Compared to the results in Sec. IIII Bl extending the fit 



Q 2 

(GeV/c) 2 


Experiment 


G% 


G s ,„ 

M 


v 2 Iv 


Prob 

(%) 


<0.11 (GcV/c) 2 


— 0.008(0.016) 


0.26(0.22) 


4.7/5 


45.5 


0.122 


G° 


— 0.007(0.016) 


0.27(0.22) 


4.9/6 


55.4 


0.128 


G° 


-0.006(0.016) 


0.30(0.22) 


8.4/7 


29.7 


0.136 


G° 


-0.006(0.016) 


0.33(0.21) 


9.9/8 


27.2 


0.144 


G° 


-0.007(0.016) 


0.30(0.21) 


13.0/9 


16.1 


0.153 


G° 


-0.007(0.016) 


0.30(0.21) 


13.0/10 


22.1 


0.164 


G" 


-0.008(0.016) 


0.29(0.21) 


13.8/11 


24.4 


0.177 


G° 


-0.010(0.016) 


0.24(0.21) 


19.5/12 


7.7 


0.192 


G° 


-0.010(0.016) 


0.22(0.21) 


20.3/13 


8.9 


0.210 


G" 


-0.012(0.016) 


0.19(0.21) 


22.1/14 


7.6 


0.230 


PVA4 


-0.013(0.016) 


0.22(0.20) 


22.3/15 


10.1 


0.232 


G° 


-0.014(0.016) 


0.21(0.20) 


23.1/16 


11.0 


0.262 


G° 


-0.012(0.016) 


0.24(0.19) 


24.7/17 


10.2 


0.299 


G" 


-0.014(0.016) 


0.23(0.19) 


26.0/18 


10.0 


0.344 


G° 


-0.014(0.016) 


0.24(0.19) 


26.1/19 


12.7 


0.410 


G° 


-0.013(0.015) 


0.27(0.19) 


27.1/20 


13.1 


0.477 


HAPPEx-H-99 


-0.015(0.015) 


0.25(0.19) 


28.6/21 


12.4 



TABLE V: Results of the first order fit (G% tx 
Q 2 ,G S M =constant). Rows are ordered byt the Q 2 of the 
data. Each row gives the results of the fit that includes all 
data in Table 11111 up to, and including, the given Q 2 . Row 
"<0.11 (GeV/c) 2 " represents the fit with low Q 2 data in- 
cluded up to PVA4-H-b. Columns "x. 
the reduced \ 2 an d the probability ( 
fit. 



2 

rain 



/v" and "Prob" are 
H P x 2d X 2 ) of the 



68.27% of confidence interval of that parameter j32|. On 
the other hand, for the two parameters [G S E and G S M ) 
that are jointly determined, the 68.27% confidence re- 
gion is instead defined by A% 2 = 2.3 contour Q. To 
demonstrate the precision of the fit, we plot the 68.27% 
(Ax 2 = 2.30) and 95% (Ax 2 = 5.99) joint confidence 
levels in Fig. Q] as the yellow (light) and gray (dark) el- 
lipses, respectively. (G S E ,G S M ) = (0,0) yields a Ax 2 of 
4.81 (/~ P x 2dx 2 = 9.0% where P % 2 is the x 2 probability 
distribution function [IH), which is depicted as the black 
cross in the figure. 



C. Inclusion of higher Q 2 data from G°, PVA4, and 
HAPPEx 

The analysis presented so far has focused on the vicin- 
ity of Q 2 =0.1 (GeV/c) 2 . It is desirable to extend the 
analysis by including data at higher Q 2 from G°, PVA4, 
and HAPPEx. As a first attempt, the previous assump- 
tion that G S E oc Q 2 and G S M is a constant was adopted. 
We shall refer to such a fit as the first order fit. Using the 
same x 2 construction as in Eqns. [TT1 and ll21 we started by 
fitting the data in TableHHwith Q 2 <0.11 (GeV/c) 2 (up 



to higher Q 2 gives consistent values for G S E and G S M . One 
also notices that for Q 2 beyond 0.164 (GeV/c) 2 , the fit 
quality (Xmin/ l/ anc ^ X 2 probability) deteriorates signifi- 
cantly, implying that our lowest order model is no longer 
able to capture the true Q 2 variation in these two quan- 
tities. 

To better characterize the data at higher Q 2 , therefore, 
one needs to introduce higher order Q 2 dependence to G S E 
and/or G S M . For y / Q 2 sufficiently below the mass of the 
kaon, chiral perturbation theory provides a systematic 
framework for characterizing the Q 2 -dependence of nu- 
cleon form factors. In this context the strange magnetic 
moment, fj, s = G S M (Q), arises at chiral order p 2 . The sub- 
leading Q 2 -dependence of G S M (Q 2 ) - sometimes called 
the strange magnetic radius - is nominally 0(p ) and is 
determined by a combination of chiral loop contributions 
and a corresponding LEC. Loop effects also generate an 
order p 3 contribution to the st rang e magnetic radius that 
is free of any new parameters [33J . However, this contri- 
bution is substantially canceled by the 0(p 4 ) loop cor- 
rections, resulting in a strong dependence on the 0(p 4 ) 
LEC [33. In contrast, G* E starts off at chiral order p 3 , 
while the sub-leading Q 2 -dependence (and corresponding 
LEC) arises at 0{p 5 ). 

Based on these considerations, we extend the previous 



7 



analysis to include all known, sub-leading Q 2 -dcpcndence 
in G S M E through 0(p i ). In practice, doing so amounts 
to including one new constant in our fit associated with 
the strange magnetic radius. Since we are interested in 
the implications for the values of G S M E at Q 2 = 0.1 
(GeV/c) 2 , we expand about this value of Q 2 rather than 
about Q 2 — 0, viz, 

G S E (Q 2 ) = G%^, G S M (Q 2 ) = G S M + M ' S (Q 2 - 0.1) , 

(14) 

where G S E and G S M again are used to represent their val- 
ues at Q 2 =0.1 (GeV/c) 2 . This will be referred to as 
the second order fit. As in the previous section, each 
row in Table Mil now leads to a new constraint in the 
form of Eqn. [Til with G S E and G S M parametrized by 
Ean.mi Again using Eqns. [TT1 and [T2l we constructed a 
proper y 2 function and performed the least square fit for 
(G E , G S M , fx' s ) by including higher Q 2 data. 

The results of this procedure are summarized in Ta- 
ble I VII Several observations can be made from the ta- 



Q 2 

(GcV/c) 2 






(GcV/c) -2 


x 2 mi J» 


Prob 

(%) 


<0.11 


-0.002(0.017) 


0.37(0.25) 


-23.7(27.5) 


3.9/4 


41.3 


0.122 


-0.006(0.017) 


0.30(0.24) 


-5.9(19.9) 


4.8/5 


43.6 


0.128 


-0.011(0.017) 


0.20(0.23) 


15.7(13.9) 


7.1/6 


30.7 


0.136 


-0.011(0.016) 


0.21(0.23) 


15.7(9.5) 


7.1/7 


41.4 


0.144 


-0.007(0.016) 


0.30(0.22) 


0.1(6.9) 


13.0/8 


11.0 


0.153 


-0.007(0.016) 


0.30(0.22) 


-0.1(5.3) 


13.0/9 


16.1 


0.164 


-0.006(0.016) 


0.32(0.22) 


-2.3(4.0) 


13.5/10 


19.9 


0.177 


-0.005(0.016) 


0.35(0.22) 


-5.8(2.7) 


15.0/11 


18.4 


0.192 


-0.006(0.016) 


0.33(0.22) 


-4.2(2.0) 


15.7/12 


20.6 


0.210 


-0.006(0.016) 


0.33(0.22) 


-3.6(1.4) 


15.9/13 


25.3 


0.230 


-0.012(0.016) 


0.35(0.21) 


-1.5(1.0) 


20.1/14 


12.6 


0.232 


-0.013(0.016) 


0.35(0.21) 


-1.6(0.9) 


20.2/15 


16.5 


0.262 


-0.012(0.016) 


0.29(0.21) 


-0.4(0.7) 


24.4/16 


8.2 


0.299 


-0.012(0.016) 


0.31(0.21) 


-0.6(0.6) 


24.8/17 


10.0 


0.344 


-0.012(0.016) 


0.28(0.20) 


-0.3(0.5) 


25.7/18 


10.8 


0.410 


-0.014(0.016) 


0.25(0.20) 


0.1(0.3) 


27.0/19 


10.4 


0.477 


-0.014(0.015) 


0.28(0.20) 


-0.1(0.2) 


28.4/20 


10.1 


<0.11 and 0.477 


-0.008(0.016) 


0.26(0.22) 


-0.3(0.3) 


4.7/5 


45.7 



TABLE VI: Results of the second order fit (G% oc Q 2 and 
G S M linear in Q 2 ). See the caption of TablelVlfor a description 
of the table content. Row "< 0.11" represents the fit using 
low Q 2 data up to PVA4-H-b. The data at Q 2 =0.230 and 
0.477 (GeV/c) 2 are from PVA4 and HAPPEx-H-99, respec- 
tively, and all others are from G°. Uncertainties of G a E , G S M , 
and fi' g are "marginalized" la uncertainties. The last row 
represents the fit by using only the low Q 2 (<0.11 (GeV/c) 2 ) 
and HAPPEx-99 data. 

ble. First, compared to the first order fits, the additional 
free parameter (i' s does not improve the fit quality sig- 
nificantly. Reasonable fits can be obtained up to Q 2 of 
0.210 (GeV/c) 2 , beyond which the "flexibility" of our fit 
model again seems to be inadequate to describe the data. 
Second, the best values of G S E and G S M are very simi- 
lar to those obtained from the first order fit (Table |V| . 



Third, the uncertainties of G S E and G S M , as compared to 
those from the first order fit, are slightly larger due to 
the additional parameter fi' s . Fourth, if we ignore the fit 
quality, and simply examine the mean and uncertainty 
of f/ s , it is large and uncertain until the fit range goes 
up to 0.262 (GeV/c) 2 . This is expected, since it is dif- 
ficult to determine the slope parameter with insufficient 
"lever arm". Also, fits beyond 0.262 (GeV/c) 2 suggests 
a gentle fx' s . (To illustrate this point, in the last row in 
Table I VII wc show the fit results by using only the low 
Q 2 (<0.11 (GcV/c) 2 ) and HAPPEx-99 data, which yield 
a consistent picture as described above.) However, the 
fit quality for Q 2 > 0.210 (GeV/c) 2 prevents us from 
making a strong statement about fi' 8 here. 

For completeness, we also investigated the impact of 
including a Q A term in G S E (corresponding to chiral or- 
der p 5 ). The resulting fit quality does not improve sub- 
stantially with the inclusion of data with Q 2 beyond 
~0.2 (GeV/c) 2 , and both second order parameters in the 
form factors are poorly constrained. Nevertheless, the re- 
sulting G S E M are consistent with those obtained above. 

Based on these considerations, we choose to use the 
first order fit up to Q 2 =0.164 (Table [VJ) as our final 
results in this analysis: 

G% = -0.008 ±0.016, G S M = 0.29 ± 0.21 , (15) 

with a correlation coefficient of —0.85 between the two, 
and Xmin/ l/ = 13.8/11. Again, the uncertainties of the 
two form factors are marginalized la uncertainties cor- 
responding to the projections of A% 2 = 1 contour. This 
result is in good agreement with Eqn. I13[ for which only 
the lowest three Q 2 bins of G° data were included. For 
G S M alone, the one-side confidence integral for a negative 
G S M ([-oo,0]) is 12.3%. 



IV. CONCLUSION 

A combined analysis of the world PV electron scat- 
tering data has been performed to extract the nuclcon 
strange electric and magnetic form factors G S E and G S M 
at Q 2 = 0.1 (GeV/c) 2 . Our treatment is similar to that 
of Ref . , but utilizes all available low Q 2 data includ- 
ing the recent HAPPEX results [13, [H| and incorporates 
the uncertainties in the electromagnetic and axial form 
factors. Wc find that the agreement among different mea- 
surements is good and we obtain fits with acceptable x 2 - 
Using a simple paramctrization of the Q 2 variation in 
G S E and G S M , a satisfactory global fit can be obtained up 
to a Q 2 ~ 0.2 (GeV/c) 2 . At Q 2 = 0.1 (GeV/c) 2 , the 
confidence integral for G S M < is 12.3%, so substantially 
negative values of G S M are highly disfavored by the fit. 
In addition, our best fit is consistent with G S E = at 
Q 2 =0.1 (GeV/c) 2 . 
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